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(Dated: February 7, 2008) 

We report the results of a search for supersymmetry (SUSY) with gauge-mediated breaking in the 
missing transverse energy distribution of inclusive diphoton events using 263 pb~^ of data collected 
by the D0 experiment at the Fermilab Tevatron Collider in 2002-2004. No excess is observed 
above the background expected from standard model processes, and lower limits on the masses of 
the lightest neutralino and chargino of about 108 and 195 GeV, respectively, are set at the 95% 
confidence level. These are the most stringent limits to date for models with gauge-mediated SUSY 
breaking with a short-lived neutralino as the next-lightest SUSY particle. 



PACS numbers: 14.80.Ly, 12.60.Jv, 13.85.Rm 



Models involving gauge-mediated supersymmetry metry breaking is achieved by introduction of new chiral 
breaking (GMSB), originally proposed in Ref. [l| have at- supermultiplets, called messengers, which couple to the 
tracted much attention i2||. In GMSB models supersym- ultimate source of supersymmetry breaking, and also to 
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the SUSY particles. The phenomenology of these mod- 
els is rich and strikingly different from that of gravity- 
mediated SUSY models. 

For GMSB models, the gravitino (with a mass less than 
~ keV) is the lightest SUSY particle (LSP), and the 
phenomenology of these models is therefore determined 
by the next-to-lightest SUSY particle (NLSP), which can 
be either a neutralino or a slepton. In the former case, 
which is considered in this paper, the NLSP decays into 
a photon and an LSP, and the signal of interest, assum- 
ing i?-parity conservation j^, is a final state with two 
photons and large missing transverse energy (J^t)- 

The model we consider is a minimal GMSB with a neu- 
tralino as the NLSP, referred to as Snowmass Slope SPS 
8 01 . This model has only one dimensioned parameter 
A that determines the effective scale of SUSY breaking. 
The minimal GMSB parameters correspond to a messen- 
ger mass Mm = 2 A, the number of messengers — 1, 
the ratio of the vacuum expectation values of the two 
Higgs fields tan/? = 15, and the sign of the Higgsino 
mass term /i > 0. The lifetime of the neutralino is not 
fixed by this model line, and is assumed to be sufficiently 
short to result in decays with prompt photons. Current 
lower limits on the GMSB neutralino mass for somewhat 
similar model parameters are 65, 75 and 100 GeV, from 
the CDF 0, D0 1^ and CERN LEP collaborations 0, 
respectively. 

We search for SUSY production in pp collisions at 
= 1.96 TeV at the Fermilab Tevatron Collider. 
The D0 detector comprises a central tracking system 
in a 2 T superconducting solenoidal magnet, a liquid- 
argon/ uranium calorimeter, and a muon spectrome- 
ter The tracking system consists of a silicon mi- 
crostrip tracker and a scintillating fiber tracker and pro- 
vides coverage for charged particles in the pseudorapidity 
range \ri\ < 3. The calorimeters are finely segmented and 
consist of a central section (CC) covering jryj < 1.1, and 
two end calorimeters (EC) extending coverage to jryj « 4, 
all housed in separate cryostats . Scintillators installed 
between the CC and EC cryostats provide sampling of 
developing showers for 1.1 < \t]\ < 1.4. The electromag- 
netic (EM) section of the calorimeter has four longitu- 
dinal layers and transverse segmentation of 0.1 x 0.1 in 
r] — 4> space (where (f) is the azimuthal angle), except in 
the third layer, corresponding to EM shower maximum, 
where it is 0.05 x 0.05. The data sample was collected 
between April 2002 and March 2004, using inclusive sin- 
gle electromagnetic (EM) and di-EM triggers. The inte- 
grated luminosity of the sample is 263 ±17 pb~^. 

Photons and electrons are identified in two steps: first, 
selection of the EM clusters, and then their separa- 
tion into photons or electrons. EM clusters are se- 
lected from calorimeter clusters by requiring that (i) at 
least 90% of the energy be deposited in the EM sec- 
tion of the calorimeter, (ii) the calorimeter isolation 
variable (/) be less than 0.15, where / = [i?(ot (0.4) — 
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FIG. 1: The l^j- distribution for the diphoton and back- 
ground samples. Also shown is the expected distribution for 
the GMSB point with A = 80 TeV, muhiplied by a factor of 
ten. 



Eem{0.2)]/Eem{Q-2), where Etotj OA) is the total shower 
energy m a cone of radius 7^ = VT^^tF+WF = 0.4, 
and Eem{0.2) is the EM energy in a cone TZ = 0.2, (iii) 
the transverse and longitudinal shower profiles be consis- 
tent with those expected for an EM shower, and (iv) the 
scalar sum of the px of all tracks originating from the pri- 
mary vertex in an annulus of 0.05 <TZ < 0.4 around the 
cluster be less than 2 GeV. The cluster is then defined as 
an electron if there is a reconstructed track pointing to 
it and a photon otherwise. Jets are reconstructed using 
the iterative, midpoint cone algorithm with a cone 
size of 0.5. is determined from the energy deposited 
in the calorimeter for I77I < 4 and is corrected for jet and 
EM energy scales. 

We select 77 candidates by requiring events to have 
two photons each with Et > 20 GeV and pseudorapidity 
1 77 1 < 1.1. To suppress events with mismeasured l^x^ 
apply the following requirements. We reject any event 
when the difference in azimuth {A(f>) between the highest 
Et jet (if jets are present) and the direction of the l^x is 
more than 2.5 radians, or if the A(j) between the direction 
of the ^^^d either photon is less than 0.5 radians. 
These selections yield 1,909 events (77 sample), out of 
which 1,800 have 1^^ < 15 GeV and two have l^x > 
40 GeV. The two events constitute the 77^7- sample. 

The main backgrounds arise from standard model pro- 
cesses with misidentified photons and/or mismeasured 
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-^•r- The background from processes with no inherent 
Iprp (multijet events, direct photon production, Z — > ee, 
etc.) is estimated using events with two EM clusters that 
satisfy photon- identification criteria (i) and (ii), but fail 
the shower-shape requirement (iii). These events, called 
the QCD sample, must pass the same trigger and other 
selections that define the 77 sample. They have char- 
acteristics similar to the background in the 77 sample 
and in particular are expected to have similar res- 
olution. This assumption was checked by varying the 
selection criteria and comparing the Ipq^ distribution in 
the QCD sample to that in Z ^ ee events. The QCD 
sample comprises 18,437 events, with 17,379 events hav- 
ing $j, < 15 GeV, and 27 events with $j, > 40 GeV. We 
estimate the background in the 77-^7^ sample resulting 
from mismeasurement of ]^rp by normalizing the number 
of QCD events to that of the 77 sample for < 15 GeV. 
This yields 2.8 ± 0.5 events with > 40 GeV, with un- 
certainty dominated by the statistics of the QCD sample. 

The other sources of background correspond to events 
with genuine in which an electron is misidentified 
as a photon, for example from W+'Y events (where 
'7' denotes both true photons and jets misidentified as 
photons), and from Z t^t^ e+e^ + X and 
tt — > e+e~ + jets production. We estimate this contri- 
bution using the e7 sample which has the same trigger, 
kinematic, and EM identification requirements as the 77 
sample. This sample contains 889 events, 782 events with 
$T < 15 GeV and 15 events with $t > 40 GeV. To esti- 
mate the contribution of such events to the 'jjI^t sample, 
we first subtract the QCD background component of the 
67 sample. This is done by normalizing the QCD sam- 
ple to the 67 sample for l^j- < 15 GeV. Then, using the 
probability for an electron to be misidentified as a pho- 
ton (measured using Z — > ee events to be 0.064 ± 0.004), 
we estimate this background to be 0.9 ± 0.2 events with 
statistically dominated uncertainty. Therefore the total 
expected background to the jjI^t sample is 3.7 ± 0.6 
events. The distributions for the 77 sample, back- 
ground without genuine ^j^, and the total background 
are shown in Fig. ^ together with an expected distribu- 
tion from the Snowmass Slope model with A = 80 TeV, 
the latter multiplied by a factor of ten for clarity. 

To estimate the expected signal, we generated Monte 
Carlo (MC) events for several points on the Snowmass 
Slope (see Table D), covering the neutralino mass range 
from 72 GeV, somewhat below the existing limits 
to 116 GeV. We used isajet 7.58 [uj to determine SUSY 
interaction eigenstate masses and couplings, pythia 
6.202 ^2 ■^^s used to generate the events after de- 
termining the sparticle masses, branching fractions and 
leading order (LO) production cross sections using the 
CTEQ5L jlj] parton distribution functions (PDF). MC 
events were processed through full detector simulation 
and reconstruction, and processed with the analysis pro- 
gram used for the data. 



The dominant contributions to the cross section 
are from production of lightest charginos (xtxi) ^'^d 
chargino-second neutralino pairs (x2xf )■ The total cross 
section in TablelHis calculated to leading order in pythia 
for GMSB SUSY production. The "X-factor" used to 
account for higher-order corrections is applied to esti- 
mate the next-to-leading-order (NLO) cross section. The 
values of the /iT-factor in the table are taken from Ref. 
[1^ . The sources of error on signal efficiency include un- 
certainty on photon identification (4% per photon), MC 
statistics (5%), and choice of PDF (5%). 

Since the observed number of events is in good agree- 
ment with that expected from the standard model, we 
conclude that there is no evidence for GMSB SUSY in 
our data. To calculate the upper limit on the production 
cross section for each sampled point on the Snowmass 
Slope, we use a Bayesian approach with a flat prior 
for the signal cross section. The calculation takes into 
account uncertainties on the expected number of back- 
ground events, efficiency, and luminosity. The selection 
> 40 GeV for the signal sample leads to the best 
expected limit, given the predicted background and ex- 
pected signal distributions. Our limits are shown in Table 
n and plotted in Fig. El together with the expected signal 
cross section. The upper limit on the cross section is be- 
low the expected value for A < 79.6 TeV, corresponding 
to lower limits on gaugino masses of m~+ > 194.9 GeV 
and m^o > 107.7 GeV. The expected limit, given the pre- 
dicted number of background events, is A > 74.5 TeV. 
We find that the gaugino mass limits depend only slightly 
on the parameters of the minimal GMSB. We have con- 
sidered models with values of tan f3 and different from 
the Snowmass Slope, and arrive at very similar results as 
detailed by Table HD 

To summarize, we searched for inclusive high-Ex 
diphoton events with large missing transverse energy. 
Such events are predicted in supersymmetric models with 
low-scale gauge-mediated supersymmetry breaking. We 
find no excess of such events, and interpret the result as 
a lower limit on gaugino masses. For a representative 
point in the parameter space, we determine that at a 
95% confidence level, the masses of the lightest chargino 
and neutralino are larger than 195 and 108 GeV, respec- 
tively. These are the most restrictive limits to date for 
the Snowmass Slope model. 
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tions, and acknowledge support from the Department of 
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Agency for Atomic Energy and RF President Grants Pro- 
gram (Russia), CAPES, CNPq, PAPER J, FAPESP and 
FUNDUNESP (Brazil), Departments of Atomic Energy 
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TABLE I; Points on the Snowmass Slope: their cross sections, efficiencies and cross-section hmits. 



A, TeV 


m^o.GeV 


m-+,GeV 


LO 1 


A'-factor 


Efficiency 


95% C.L. Limit, pb 


55 


71.8 


126.3 


0.735 


1.236 


0.092 ± 0.009 


0.184 


60 


79.1 


140.2 


0.468 


1.227 


0.100 ±0.009 


0.170 


65 


86.4 


154.3 


0.301 


1.217 


0.111 ±0.011 


0.153 


70 


93.7 


168.2 


0.204 


1.207 


0.124 ±0.012 


0.137 


75 


101.0 


182.3 


0.138 


1.197 


0.137 ±0.013 


0.124 


80 


108.2 


196.0 


0.094 


1.187 


0.149 ±0.014 


0.114 


85 


115.5 


209.9 


0.066 


1.177 


0.154 ±0.015 


0.110 



TABLE IL Limits on the Snowmass Slope and two other 
GMSB models. 

Fixed parameters 95% CL lower limits 

Mm/ A tan/3 TVs sign(^) A m^o m-+ m^o 

2 15 1 + 706 IOTTt 1949 19^9" 

2 5 1 + 79.5 106.0 191.6 193.3 

10 5 2 + 44.0 111.4 196.0 198.7 
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FIG. 2: Predicted cross sections for the Snowmass Slope 
model vs A in leading order (thin solid line with crosses), 
multiplied by the iC'-factor (thin dashed line), and the 95% 
C.L. limits (solid line). 
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